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Abstract. We attempt to employ measures of average fidelity and fidelity deviation 
■ in characterizing a quantum operation when realizing the universal NOT gate for a 

^ ■ single qubit. Based on these measures, we evaluate quantum operations with respect to 

optimality and universality. We show that one-qubit operations have a sharp trade-off 
relation, i.e. a linear relation between average fidelity and fidelity deviation. Two- 
qubit operations for the same task show a looser trade-off relation between the two 
quantities. We show that the genuine universality comes to hold for operations of 
more than 2 qubits, while those of even more qubits are nevertheless beneficial to 
compensating imperfection of control. In addition, we take into account operational 
noises which contaminate quantum operation in realistic circumstances. We show that 
the operation recovers from the contamination by a feedback procedure of differential 
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I evolution. It is discussed that our feedback scheme is also applicable to finding an 

optimal and universal operation of NOT. 



X: 



PACS numbers: 03.65.-w, 03.67-a, 03.67.Lx, 03.67.Pp 



2 



1. introduction 

Quantum information offers advantages in variety tasks over classical counterparts, by 
virtue of fundamental properties of quantum physics pLj. Quantum theory imposes, on 
the other hand, certain restrictions on quantum tasks [2]. For example, an arbitrary 
quantum state cannot be cloned, called no-cloning theorem, so that the superluminal 
(i.e. faster than light) communication via entanglement is prohibited [3l 0]. Another 
quantum task of universal-NOT (U-NOT) that transforms an arbitrary input state to its 
orthogonal is also restricted by quantum theory, while its classical task NOT is perfectly 
realized by bit-flip operation [5l E] . This is because U-NOT cannot be implemented by 
a unitary operation but by an anti-unitary operation, violating the conditions of trace 
preservation and complete positivity that a physical procedure obeys [3 [HI [9] . 

An approximate realization to quantum task U-NOT can nevertheless be done by 
a physical operation assisted by ancillary qubits |5l |6]. The approximate operation 
is the most optimal when it yields the average fidelity 2/3 [5]. Such an optimal 
operation has extensively been studied for the last decade both theoretically [5l|6l[T0] and 
experimentally [HI [121 [H] in order to clarify capabilities and limitations of quantum 
information processing. In another perspective, the optimal operation of U-NOT is 
closely related to other important quantum tasks such as quantum cloning, quantum 
state estimation, and entanglement test [5l [101 HH HH US]. In particular, U-NOT is 
equivalent to the transposition by some unitary transformation [151 [13] • This implies 
that the optimal operation of U-NOT enables to approximately test if a quantum state 
is entangled with negative partial transposition [161 E] ■ 

An approximate physical operation has been evaluated in terms of the fidelity 
between its output state and the target of the task. In particular, the average fidelity 
over all possible input states has been employed as an optimality measure of the 
operation to the task. However the average fidelity itself tells nothing about universality 
of operation, the condition that the task is performed equally for all possible input states. 
In a theoretical side, universality can be imposed on a quantum operation by requiring 
it to result in an equal fidelity for all input states. On the other hand, such requirement 
is nontrivial in experiments where imperfections of control and noises by environment 
arise. It is thus desired to consider a measure to quantify the state-dependence of 
operation over input states. As such a measure, we employ fidelity deviation which is 
defined by the standard deviation of fidelity over all possible input states. 

In this paper, we attempt to employ fidelity deviation A and average fidelity F 
in characterizing a quantum operation of U-NOT and to represent the operation by 
a point on two-dimensional space of F and A. By this way we analyze a quantum 
operation in terms of optimality and state-dependence. We show that interestingly 
some one-qubit operations can already reach the average fidelity of 2/3, the maximum 
by three-qubit U-NOT operation [5], whereas they exhibit strong state-dependence. It 
is proved that a quantum operation, no matter how many qubits are involved in, can 
not exceed 2/3 of average fidelity. We show that one-qubit operations have a sharp 
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trade-off relation, i.e. a linear relation on the space of (F, A). Two-qubit operations 
have a less sharp trade-off relation, including the one-qubit relation as an upper bound. 
They can be the most "universal" with A = 1/3-^/5 ~ 0.15. The genuine universality of 
A = holds for n-qubit operations with (n — 1) ancillary qubits if n > 3. We discuss 
that operations of more than 3 qubits are beneficial against imperfection of control. 
In addition, considering some realistic circumstances, we take into account operational 
noises which contaminate quantum operation once optimized. We find a case that such a 
polluted operation is far from the universality no matter how close its average fidelity is 
to the maximum of 2/3. To protect an operation against operational noises, we suggest 
a feedback scheme of using a differential evolution, showing that our scheme recovers 
the operation from the contamination as far as the noises fluctuate slowly compared 
to the operation. It is discussed that our scheme of feedback is applicable to find an 
optimal operation of U-NOT with no a priori knowledge except the number of qubits. 

This paper is organized as follows. In Sec. II, we introduce our approach for the 
optimality and the universality by employing average fidelity and fidelity deviation. 
Sec. Ill is devoted to analyses of quantum operations for task U-NOT on the two 
dimensional space of average fidelity and fidelity deviation. In Sec. IV, we investigate 
effects of operational noises on optimal operation of U-NOT and we suggest a feedback 
scheme to cure the contaminated operation by noise. Remarks are given in Sec. V. 

2. Average fidelity and fidelity deviation 

A task is realized by a physical operation O that transforms an input state to its 
target state Some tasks can not be ideally realized as restricted by quantum laws, 
such as quantum cloning [3]. It is thus desirable to find an approximate but optimal 
operation as close as possible to a given task. 

The quality of a found operation O is commonly quantified by a quantum fidelity 
/, that is defined by the transition probability between the output state O and the 
target state of an input state 

/[Vl/] = |(V1>,|6|V]>)|2. (1) 

For a given operation, the fidelity varies in general on input states. If not, the 
quantum operation is said to be universal. In that sense, by the universal operation, 
the task can be performed equally for all possible input states. However, in some 
realistic circumstances, it is difficult to achieve the universality, due to noise (s) during 
the physical process and/or an unavoidable interaction with environment. It is thus 
necessary to introduce a measure of state-dependence to quantify how much the task 
depends on the input states and to examine whether to reduce the state-dependence 
by altering experimental parameters. In this section, we propose to employ the average 
fidelity for quantifying the optimality of an operation and the fidelity deviation for 
quantifying the state-dependence. 
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The average fidelity F is defined as 



(2) 



where the integral is over all possible input state |\1') and d'^ is a normalized Haar 
measure, J d'^ = 1. The measure F quantifies on average how well operation O 
transforms input state to their target states {"^t)', the value F = 1 implies the 
task is perfectly performed for all possible inputs, while F = 1/2 does a random task. 
We quantify the state-dependence of the operation O in terms of the standard deviation 
of/, 

nl/2 

We call A a fidelity deviation. The fidelity deviation A has the minimum of if / = F 
for all input states and otherwise it increases. Note 

1 



A^ < / d^fl"^] -F' = F{1 -F) < 



4' 



(4) 



where the last equality holds when F = 1/2. Thus, A is bounded as < A < 1/2. 
By the two measures F and A, we characterize a task operation, as a point on the 
two-dimensional space of {F,A). 



3. Characterization of U-NOT operations on the space of (F, A) 

In this section, we consider approximate operations for implementing task U-NOT and 
place them on the space (F, A). An optimal operation of U-NOT was found among 
three-qubit operations [51 E] . It is questioned whether there exist any operations of U- 
NOT among one- or two-qubit operations. We try to answer this question and generalize 
to arbitrary number of qubits. 



3.1. One-quhit operations for U-NOT 

An input state of all possible pure states is given in the Bloch representation by 

A„ = |vl/)(vl/| = l(i + aV), (5) 

where 1 is the identity operator, a = {ax, ay, az)^ = tr (cr't'pin) is a Bloch vector of unit 
norm in 3-dimentional real vector space M^, and cr = (a^;, a^,, a^)^ is a vector with its 
components being Pauli operators aj (j = x, y, z). Note that all pure states are located 
on the surface of the Bloch sphere with |ap = 1. The task U-NOT is supposed to 
transform each input to its orthogonal state or the Bloch vector a to its antipodal 
—a: 

A4=|^^)(^1=^(i-aV), (6) 
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The state is the target state of task U-NOT. To find a physically realizable 
(approximate) operation for U-NOT, we consider an arbitrary one-qubit unitary 
operation, given by 

f ^ T \ ^ ' ^ / T N /«N 

U = exp \—i—n a j = cos — 1 — i sin (' ) 

where n = [rix, Uy, UzY" is a unit vector. The operation transforms pin to the output 
state, 

where b = Ra and R is a rotation matrix on M^. The operation U can be understood 
as a rotation R, on the Bloch vector a, of the angle d along axis n [H [18]. Note that 
the output state is also pure, i.e. |bp = 1. 

The fidelity / = tr (p^ Pout) between the output state pout and the target state p-^ 
is given by [19], 

/[a] = i(l-a^Ra). (9) 

The average fidelity over all possible input states or all Bloch vectors a on the Bloch 
surface is given by 

FiQ = j c/a/[a] = j rfa^(l-a^Ra), (10) 

where cia is the (normalized) Haar measure over the surface of the Bloch sphere |20j . 
The subscript "IQ" stands for one qubit. 

Equation ( !T0|) is evaluated in a spherical coordinate system, where a"^ = 
(sin ^ cos 0, sin 6^ sin 0, cos 6*) and (ia = sin 9d9d(j). The diagonal components aiRaai 
are integrated to be |-Rm, while the non-diagonal aiRijaj are to vanish. Alternatively, 
one may utilize Schur's lemma (Sec. 2.2 in reference [21]). 

[0,XOj]^ = ^tr(X)I,, (11) 

where 1^ is an identity matrix in ci-dimensional real vector space M"', Og is an irreducible 
orthogonal representation of an element g in a given group G, and [Fg]G denotes the 
average of Fg over all elements g E G: [Fg]G = J dgFg, where dg is the (normalized) 
Haar measure such that J dg = 1. This holds for every matrix X on M'^. By applying the 
lemma f[TTl) to the group 0(3) of 3-dimensional rotations, the second term in equation 
f[TU]) results in 

c/aa'^Ra= -tr(R), (12) 
3 

where we used the fact that every Bloch vector a is given by some rotation R from a 
certain reference z, a = Rz, and the average over the Bloch sphere is equal to that over 
the rotation group 0(3), / (iaRj|-j^^RRg(a) = / dgH^URg = |tr (R)l3. Both methods 
result in 

FiQ = ^- Jtr(R), (13) 
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where tr (R) = 2cos'i9 + 1 (see Appendix A). The maximum of Fiq is given to be 2/3 
when d = TT and the minimum is when -(9 = or 27r. It is remarkable that the maximal 
average fidehty of one-qubit operation is already equal to that of three-qubut operation 
for U-NOT O [6]. In the case, the found optimal operation is in the form of equation 
(I7j) with = 71 and n being an arbitrary unit vector. 

We investigate the state-dependence of one-qubit operations for task U-NOT. The 
square of the fidelity deviation is 



(ia / [a 



2 _ f2 



dsi (a^Ra)^ - itr(R)^ 



(14) 



To evaluate A^g in equation (fT^ . we use a generalized identity of Schur's lemma in 
equation ( ITTll to the tensor product of the two real vector spaces M'^ (g) ] 
matrix X on M"' (g) W^, 

[(O, ® Og) X (Oj ® Oj)]^ = +/3-D + 7P, 



i'^, given for each 
(15) 



where 



a 



/3 
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{d + 1) tr (X) - tr (XD) - tr (XP) 
d{d-l){d + 2) ' 

- tr (X) + {d+l) tr (XD) - tr (XP) 

d{d-l){d + 2) '' 

- tr (X) - tr (XD) + {d + 1) tr (XP) 



d{d- l){d + 2) 



Here, P is a swap matrix P (xj ® x^) = x^ ® Xj, or equivalently. 



d-l 



i,j=0 



and 



D 



J=0 




where {xj} is an orthonormal basis set in W^. Then, using the identity of equation ([T5 
we rewrite the first term in equation (fHIl . 



j da (a^Ra)^ = j da (a (g) a)^ (R (g) R) (a (g) a) 



1 
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tr (R ® R) + tr (R (g RD) + tr (R (g RPl 



(16) 



where we used the similar reasoning below equation f|T2l) . Note that tr (R (g R) 



tr(R)^ tr(R®RD) = tr(RR^) = tr (I, 



equation (]T41) is rewritten as 
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[tr (R) 



3 + tr (R' 



3, and tr (R® RP) 
1 



tr(R2). Then 
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■tr (R) 



IQ' 



(17) 
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Average fidelity F 

Figure 1. Accessible region of quantum operations for task U-NOT in terms of average 
fidelity F and fidelity deviation A. A one-qubit operation lies at a point on line OPi, 
where O is the origin. An operation assisted by one quit (or two qubits) lies inside 
or on triangle OP1P2 (or OP1P3). The (blue) dashed line stands for a mathematical 
boundary of quantum operations [see equation Q]. 



where we used the relation tr (R)^ — tr (R^) = 2 tr (R) in Appendix A. The final form 
of the fidelity deviation Aiq is given by 

AiQ = i=FiQ. (18) 

We note that this relation holds for arbitrary one-qubit operations as well as the optimal 
operations. This relation is represented by a segment OPi in the space (F, A), as 
shown in figure [H equation (fTS!) clearly shows the sharp trade-off relation between the 
optimality and the state-dependence of one-qubit operations for task U-NOT: The larger 
the average fidelity (optimality), the larger the fidelity deviation (state-dependence). 
Thus, there is no one-qubit operation that satisfies both of universality and optimality 
as the condition Aig = demands Fiq = 0, i.e. an identity operation, even though the 
maximal average fidelity is equal to that of three-qubit U-NOT. 

3.2. n-qubit operations assisted by (n — 1) qubit(s) 

We generalize one-qubit to ra-qubit operations, by employing a specific type of logic 
circuits, as seen in figure |2j In the circuit, the first qubit is the system and the 
rest of (n — 1) qubit (s) are ancillary. The case of n = 1 was investigated in the 
previous sub-section. The circuit operation consists of local unitary Vj and conditional 
unitary Uj. The local unitary operator Vj on the ancillary qubit j is defined such that 
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Figure 2. A quantum circuit for task U-NOT with (n — 1) ancillary qubits. The part 
inside the dashed small (large) box corresponds to one(two)-qubit assisting circuit. 



~ + — Vj where Vj is a real number, satisfying < f j < 1. The 

conditional unitary operator Uj acts on the systemm conditioned that ancillary qubits 
k are in the state 11);^. for all k < j. 

When the circuit operation is applied on an input state |\E') of the system and the 
states |0) of the ancillary qubits, the output state pout becomes, partially tracing over 
the ancillary qubits, 

rt— 1 ^ 



■)• 



(19) 



A:=0 



Here, positive Wk are functions of v/s, 



Wk 



(fc = 0) 
(1 < A; < n 
{k = n-l) 



2) 



(20) 



satisfying XlLo 

Wk = 



(21) 

IbP < 1, due to 



(1 - Vk)Vk+l 

(1 - t;„_2)(l - v„ 
= 1 and the unitary operators Wk are given by 
UkUk-i ■ ■ ■ Uq. 

Note that the output state pout is not necessarily a pure state, i.e. ^ ^ _ 
the entanglement created during the process. Nevertheless, it is remarkable that the 
circuit can be understood as a stochastic unitary map, in equation f|T9l) . characterized 
by the set of the local unitary operators Wk and the probability weights Wk- In other 
words, the circuit can be replaced by a stochastic circuit on the system that the unitary 
operation Wk is applied in the probability Wk- In the sense, the entanglement presented 
in the circuit is not necessarily demanded. The stochastic representation of operation in 
equation f|T9|) reduces significantly the calculations in the characterization of the average 
fidehty and the fidelity deviation. 

Consider the average fidelity F„q (where the subscript "nQ" stands for n-qubit). 
The fidelity of the output state in equation ( !T9|) to the target state is given by 



n-l 



/[a] = tr (p^pout) = ^Wfc/4a], 



(22) 



A:=0 
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where /^[a] = (l — a-^R^a) /2 and the rotation matrices Rfc 6 MP are associated with 
Wk- Each Wk and thus Rfc is given by the rotation angle -^k and axis n^., as in equation 
(^^. The average fidehty F„q is given as 

// n—1 \ n—1 

da I ^Wfe/fc[a] J = ^WkFiQ^k, (23) 
\A:=0 / k=0 

where -Fiq ^ is the average fidehty by a one-qubit operation Wk, as in equation ffTOj) . It 
is clear that < F^q < 2/3, as F„q is a statistical mean of average fidelities of one-qubit 
operations. The maximum F^q = 2/3 is attained when Fiq^^ are all equal to 2/3, or 
equivalently -^k = tt for all k. This result holds for an arbitrary number of ancillary 
qubits. It thus seems that increasing the number of ancillary qubits does not improve 
the average fidelity or optimality for task U-NOT. 

The square of the fidelity deviation A^g is given from equations ( l22i) and ( 123|) as 

//n-l \ 2 n-l 

rfa J2 ^'^/^[^l - ^nQ = E ^f'^iCf'i^ (24) 

\fc=0 / k,l=0 

where Cki are elements of covariance marix C, defined by 

Cki = J cia/fc[a]/Ja] - F,Q,kFiQ,i. (25) 

Note that C is symmetric, i.e. Cki = Cik- Each element of C is bounded, as shown in 
Appendix B, by 

Ckk = ^lQ,k 

(26) 

— 2^^Q,k^lQ,l ^ < ^lQ,k^lQ,l, 

where Aiq^k is the fidelity deviation of one-qubit operation Wk- The equality for the 
lower bound holds when the two rotation axes and are orthogonal to each other, 
i.e. n^ri; = 0, and the upper bound is reached when n^. and are parallel or anti- 
parallel, i.e. n^n/ = ±1. By equation (l26l) . the fidelity deviation A^g in equation ( l24l) 
is upper bounded. 



'n-l ^ 2 



^Iq< E^^^IQ'O =l^nQ, (27) 
\A:=0 / 

where we used equations f|T8l) and fl23l) . The equality holds when n^^n^ = ±1 for all 
pairs oi k ^ I. The lower bound of A^g is given as 

n—1 ^ n—1 2 n—1 

^Iq > ^I^Iqm - 2 5Z WkWi/\^Q,k/^^Q,i > (28) 

fc=0 fc^i fc=0 

where we used equation fl26l) and the inequality, Xlfc^^i ("u^fc^iQ.fc — wiI\\Q^if > 0. The 
two equalities successively hold when n^n/ = and WkAiq^k = const for all pairs of 
k^l. 
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Assisted by single ancillary qubit. - Based on the above results, let us consider 
two-qubit operations. The circuit is depicted inside the small dashed box in figure [2j 
The stochastic probabilities Wk {k = 0,1) are given by equation ( l20l) . wq = vi and 
Wi = 1 — ^1. The average fidelity F2Q ranges from to 2/3. When n^ni = and 
woAiQfi = wiAiQ^i, the lower bound of A^g in equation ([281), (1/2) X]fc=o,i ^fc^iQ.fc' 
is attained and it is equal to {i/4:){^f,^Q^WkAiQ^kY = -^2^q/20. Thus, the following 
inequalities hold, 

^ F2Q < A2Q < ^F2Q. (29) 



2^5 - V5 

This implies a trade-off relation of F2Q and A2Q for two-qubit operations, as represented 
by the triangle OP1P2 in figure [TJ The trade-off relation in equation f l2I?]) is looser 
than one-qubit operations in the sense that for given average fidelity F we can always 
find an operation of two qubits which is less state-dependent than that of one qubit. 
The most optimal and "universal" operation is given when the operation Wq and Wi 
satisfy ■})() = i^i = ir for their angles and n^ni = for their axes, and the stochastic 
probabilities wq = wi = 1/2. The state-dependence is reduced to A2Q = 1/3 a/5 ~ 0.15 
for the optimal operations of F2Q = 2/3. We note that the circuit operations we have 
considered include all possible two-qubit operations and the current results hold in 
general as far as two qubits are involved. 

Assisted by two ancillary qubits. - Consider three-qubit operations, as shown in the 
large dashed box in figure El The local one-qubit unitary V2 and controlled-controlled-t/2 
operators are additionally employed for the task, and the stochastic probabilities are 
given by wo = vi, wi = (1 — Vi)v2, and W2 = (1 — fi)(l — ^2)- The lower bound A3Q = 
is reached when the three vectors n^, are mutually orthogonal and w^Aiq ^ = const, Wk. 
Thus, we arrive at the trade-off relation, for three-qubit operations, 

< A3Q < -^FsQ. (30) 

This relation is represented by the triangle OP1P3 in figure [H The most optimal and 
genuinely universal operation of F^q = 2/3 and A3Q = is attained when the stochastic 
unitary operations Wk are given by their rotation angles 'dk = and their axes mutually 
orthogonal n^^n^ = with Wk = 1/3 for all pairs of k ^ I. In terms of a stochastic map, 
the most optimal operation of U-NOT leads 

1 2 1 

Pin Pout = - {o-xPinO-x + CTyPinC^y + (T^PinO-^) = -p^^ + gPin- (31) 

This map is equivalent to the one found in references P, 122] ■ 

The result in equation (l30l) still holds for more than 2 ancillary qubits. Our analyses 
show that it is important to employ both indicators of the average fidelity and the 
fidelity deviation to evaluate a quantum operation of task U-NOT, because there exist 
operations whose average fidelity F are close to | but fidelity deviation A may be 
arbitrarily large, as implied by the line P1P3 of figure [H It is understood that such 
situation could be a case in experiments, as in references [T2|1T3]. which will be discussed 
further in Sec. IV. 
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Before closing this section, remind that the three-qubit operations we have 
considered are the specific, as in the circuit, figure [2l One might question if there exist 
any three-qubit operation whose average fidehty is larger than 2/3 when sacrificing the 
universality. This question is worth to investigate as the universality was assumed in 
the previous works [HI El [121 [13] • However, this is not the case. Consider an arbitrary 
three-qubit operation, 

/ Mo,o Mo,l ■ ■ ■ M0,7\ 
Ulfi ■ ■ ■ Ml,7 



arb 



(32) 



\M7,0 "7,1 ■■■ ^7,7/ 

The average fidelity F^q is a function of the matrix elements Uj^k {j and A; = 0, 1, ■ ■ ■ ,7), 

^"^/' 2 2 2 2'\ 

-^SQ = 2 - g |Mo,0 + "mI +\Ui,0 + U5a\ + 1^2,0 + %,4| + |«3,0 + "7,4! j- (33) 

The unitary condition, f/arbf^irb — 1; leads to < F^q < 2/3. This proof can 
straightforwardly be generalized to arbitrary n qubit operations. 

There arises another question: Is there any advantage in using more than 2 ancillary 
qubits? The answer is affirmative: Added qubits can be used to compensate or to absorb 
imperfection of operations if any. To see this, suppose that a stochastic operation 
Wk satisfies dk = and Wk = 1/3 for all k = 0,1,2 and two rotation axes are 
not perfectly orthogonal with the angle f — « for small a, say n^ni = n^n2 = 
and njn2 = cos(| — a) ~ a, neglecting higher order terms than a"^. Then, even 
though the optimality is achieved with F^q = 2/3, the universality is broken as 
A3Q ^ 2a/3vT5 7^ 0. In such circumstance, universahty can be cured by extending 
the circuit from three to four qubits with chosen such that i)^ = tt and is at the 
opposite direction to 112 by a on the plane ni-n2, that is, n^iiQ = 0, 113 ni = n^iii ~ a, 
and 113 n2 = cos 2a 1 — I0? . By choosing the stochastic probabilities wq = w\ = 1/3, 
and W2 = W3 = 1/6, then, the state-dependence becomes to A4Q ~ up to a^, 
while keeping F4Q ~ F^q = 2/3. This example opens a possibility of recovering 
the universality without sacrificing any optimality when operations suffer from the 
imperfection. 



4. Feedback scheme to stabilize a quantum operation 

Implementing a quantum operation suffers from noise in realistic circumstance. To 
protect, we consider a feedback procedure with a differential evolution method, which 
is known as an efficient heuristic method for global optimization [23]. The adotion of 
such a feedback procedure is also beneficial when to find a quantum operation of itself. 
In this section, we introduce the differential evolution briefly and apply to the problem 
of flnding an optimal operation of U-NOT among three-qubit operations. We show that 
the feedback scheme works so well that it consistently flnds optimal operations of U- 
NOT, equivalent to the one in reference [5]. By introducing an operational noise which 
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Figure 3. Average fidelity F (left) and state dependence A (right) of the optimal U- 
NOT operation, when contaminated by operational noise of degree rj. For each 77, we 
perform 1000 simulations of Monte Carlo, averaging F and A by the sample. Error bars 
are their standard deviations. Both of F and A are degraded by the operational noise 
and such behaviors become conspicuous as increasing 77. Dashed lines are of a random 
operation, iv = 1/2 and A^ ~ 0.150, and dotted line in left graph is = 2/3 ~ 0.666. 



alters operational parameters unexpectedly, we show that the contaminated operation 
is cured by the feedback scheme as long as the noise fluctuates slowly. 

4-1- Effects of operational noise 

Unitary operations on rf- dimensional Hilbert space are parameterized by {(P' — 1)- 
dimensional real vectors p = (pi,P2, ■ ' ' iVd^-i)^ 

C/(p)=exp(-2p^G), (34) 

where G = {gi, g2, - ■ ■ , QdP-iY ^ vector whose components are SU((i) group generators 
(jj (j = 1, 2, ■ ■ ■ (i^ — 1) [211 EHl |26]. Components of p are control parameters. Such a 
unitary operation can be realized by multiport beam splitters for an optical system [27] 
or pulse sequences for nuclear magnetic resonance system [28]. Based on the analysis 
of the previous section for U-NOT, we consider three-qubit unitary operations t/(p) 
on 8-dimensional Hilbert space with 8^ — 1 = 63 control parameters. Note that the 
number of control parameters can be reduced if any restriction on quantum operations 
are imposed, even though we assume no restrictions throughout this paper. 

In the presence of noise, an operation f/(p) turns to be imperfect with fluctuation 
of p [29]. We choose a noise model in which fluctuation arises when dialing the control 
parameters pj such that 

p^p + r^e, (35) 

where e is a random stochastic error vector whose components are random between 
— TT and vr. The factor rj, normalized in [0,1], stands for the degree of inaccuracy in 
control. This type of noise is supposed to occur in implementing [/(p) and it is called 
an operational noise. 
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For U-NOT operations, we present the effects of operational noise on average fidelity 
F and fidelity deviation A in figure |3l The average fidelity F decreases and the fidelity 
deviation A increases as the degree of noise r] increases. That is, the performance 
of the operation is degraded, as expected. It is remarkable that for a small noise the 
average fidelity F remains close to its maximum 2/3 but the fidelity deviation A becomes 
rather large toward that of random operation, A^ = l/3\/5. For instance, when 77 = 0.1, 
averaged over 1000 samples, F is 0.633±0.018 which is about the value in an experiment 

[13], whereas A is rather high of 0.095 ± 0.027, compared to A,. ^ 0.150 (see figure 
|3]). In other words, F is degraded by about 25% from its maximum 3/2 to that of 
random operation, F^ = 1/2, whereas A is increased by about 65% to A^. The results 
support again the importance of the state-dependence in experimentally implementing 
a universal operation. 



4-2. Recovery from the contamination 

Our differential evolution scheme of feedback is to find a set of values for control 
parameters p for an optimal and universal operation of NOT. The differential evolution 
algorithm follows [23]. To begin with, we account A^pop operations by which we are 
to develop approximate solutions. Then, we have iVpop parameter vectors pj [i = 
1, 2, ■ ■ ■ , A^pop), each of which consists of 63 components pj^i G [— tt, vr] (j = 1, 2, ■ ■ ■ , 63). 
All these 63 x A^pop parameters are chosen initially at random. [S.l] We generate A'pop 
mutant vectors Ui according to 

h'i = Pa + D (pt - Pfc) , (36) 

where we randomly selected a, b, and c among A^pop parameter vectors as far as 
they are mutually different. The free parameter D, called a differential weight, is 
a real and constant number we choose. [S.2] After that, the parameter vectors 
Pi = {pi,i,P2,i,- ■ ■ ,P63,iV are reformed to trial vectors = {ri^i,T2,i, ■ ■ ■ ,^63,*)^ by 
the following rule: For each j, 



Tj,j i^j^i otherwise. 



(37) 



where Vj G [0, 1] is a randomly generated number and the crossover rate CR is another 
free parameter we choose in [0, 1]. [S.3] Lastly, the trial vector Tj is taken to be p, for 
the next iteration if it yields a larger fitness value than pj, and otherwise pi is retained. 
Here the fitness ^ of a given operation U{p) is defined by 

e = F - A, (38) 

where F and A are the average fidelity and fidelity deviation for U{p), respectively. 
It tells us how fit U{p) is to an optimal and universal operation of NOT [||. The steps 
[S.1]-[S.3] are repeated until the maximum iterations. 

I One might take another function of fitness. For instance, we tried ^ — F{1 — A) but the overall 
tendency did not change much. This example does not exclude existence of a certain fitness function 
which would lead to very different behaviors. 
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Figure 4. Searching an optimal and universal operation of NOT by our feedback 
scheme of differential evolution in terms of average fidelity F (left) and state 
dependence A (right). For each iteration, we perform 1000 Monte-Carlo simulations 
in averaging F and A with their error bars. As iterating, F and A steadily approach 
to their ideal optima, F = 2/3 (dashed line) and A = 0. We obtain F = 0.663 ± 0.002 
and A = 0.006 ± 0.002 in 1000 iterations. 



We perform Monte-Carlo simulations. In the simulation, we take A^pop = 10, and 
the free parameters D = 0.1 and CR = 0.03 which optimize our simulation. At every 
iteration, the fitnesses of all the operations are evaluated to select suitable parameters 
Pj for the next iteration, as described in [S.3]. We terminate the feedback procedure on 
1000 iterations. Figure H] presents the average fidelity F and the fidelity deviation A of 
the best among A^pop operations at every 20 iterations. Both F and A are statistically 
averaged by 1000 trials of simulations. As seen in figure HI F converges to its ideal 
maximum 2/3 ^ 0.667 and A also converges to zero. We obtain F = 0.663 ± 0.002 
and A = 0.006 ± 0.002 in 1000 iterations. This result shows that our feedback scheme 
can be used to search the optimal and universal (or the least state-dependent, if not) 
operation of NOT with no a priori knowledge on it, once the number of qubits is fixed 

m- 

We also perform a Monte-Carlo simulation to test if our feedback scheme is able to 
recover the operation once contaminated by the operational noise. We assume that the 
noise fluctuates slowly compared to the operation, which is the case in most experiments 
for quantum tasks [311 132]. Two cases are investigated that the abrupt fluctuation of 
noise occurs at every 50 or 100 iterations. Accounting the large degree of noise in 
equation (135|1 . we take rj = 0.5 for figure [51 Here, the operation initially optimized is 
polluted by the noise at every 50 or 100 iterations, on which the average fidelity F and 
the fidelity deviation A suddenly deteriorate close to those of random operation. As the 
feedback goes on, however, they steadily recover to those of the optimal and universal 
operation. Note that the degree of recovery depends on both frequency and degree of 
noises, as one may expect. 
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Figure 5. Recovery of contaminated operation by the feedback scheme in terms of 
average fidehty F (left) and state dependence A (right). The operational noise is 

assumed to occur at every 50 (dashed line) or 100 (solid line) iterations with the noise 
degree rj = 0.5. The dotted line in left is F = 2/3 ~ 0.666. Whenever the noise 
occurs, F and A suddenly deteriorate close to those of random operation. As iterated, 
however, they steadily recover to those of the optimal and universal operation of NOT. 
The degree of recovery depends on the frequency of noise. 




5. Remctrks 

Wc proposed to employ fidelity deviation as well as average fidelity in characterizing 
a quantum operation which is supposed to perform a quantum task of universal NOT. 
The average fidelity F represents the optimality of the quantum operation on average, 
while the fidehty deviation does the state dependence (reciprocally, "universality") of 
the quantum operation. An example is easily found for which the quantum operation 
is rather strongly state-dependent even though its average fidelity is very close to its 
maximum: The quantum operation is far from the universal. We suggested a quantum 
operation to be represented by a point on two-dimensional space of (F, A), by which 
way we analyzed the operation with respect to optimality and state-dependence. 

We showed that interestingly some one-qubit operations can reach the average 
fidehty of 2/3, the maximum by three-qubit optimal operation of task U-NOT, whereas 
they exhibit strong state-dependence. It was proved that there exists a quantum 
operation for arbitrary number of qubits such that it leads to the average fidelity of 
(but not larger than) 2/3. The one-qubit operations show a sharp trade-off relation, 
i.e. a linear relation between the average fidelity and the state-dependence. Similar 
behaviors were observed for two-qubit operations with one qubit being ancillary. The 
two-bit operations have a less sharp trade-off relation, i.e. a triangular region on the 
space of (F, A), which includes the one-qubit relation of trade-off as an upper bound. 
They can have the most "universality" of A = 1/3 -\/5 0.15. The genuine universality 
of A = was shown to hold for n-qubit operations with (n — 1) ancillary qubits as far 
as n > 3. Even though 3-qubit operations suffice to optimally perform the task U-NOT, 
it was shown that more-qubit operations are beneficial against certain imperfections 
involved in. 
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In some realistic circumstances, operational noises may arise in imperfect control of 
operation. The noises contaminate quantum operations once optimized. We emphasized 
the existence of case that such a polluted operation is far from the universality no 
matter how close its average fidelity is to the maximum of 2/3. This result supports 
again our proposal of employing the measure of state-dependence. In order to protect 
an operation against operational noises, we proposed a feedback scheme of using a 
differential evolution. It was shown that our scheme recovers the operation from the 
contamination, as far as the noises fluctuate slowly compared to the operation. We 
showed that our scheme of feedback is also applicable to find an optimal and universal 
operation of NOT with no a priori knowledge except the number of qubits. 

We expect that our proposal of employing the measures of state-dependence and 
average fidelity will be applied to other universal quantum tasks such as cloning, 
teleportation, and inseparability test. Its modifications are eligible for partially universal 
tasks which involve a subset of states. 
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Appendix A. tr (R) and tr (R^ 
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We evaluate the traces of rotation matrices on three-dimensional real vector space M^. 
For the purpose, it is useful to represent a rotation matrix R in Rodrigues' form [33] 

R = l3-sin^9S + (l-cos^9)S^ (A.l) 

where is the rotation angle and I3 is the identity matrix on M^. Here, S is the skew- 
symmetric matrix of cross product of the rotation axis n = {nx,ny,nz)'^ , defined as 

A 

(A2) 

(A.3) 
(A.4) 
(A.5) 

which depends on the rotation angle "i? but not the rotation axis n. An alternative way 
to obtain Tr(R) is to find and sum eigenvalues of R. As R has eigenvalues of {1, e^"^}, 
their summation is equal to equation ( ]A.5I) . 



where eijk is Levi-Civita symbol. The squre of S in equation (1A.2|) is written as 

= n (g) n"^ - I3. 
From equations (1A.2I) and (I A. 3 II we obtain 



tr(l3) = 3, tr(S) = 0, and tr (S^ 
Thus, the trace of R is given as 

tr (R) = 2cos?9 + 1, 



-2. 
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We now prove the relation, 

tr(R)^ -tr(R2) = 2 tr (R), (A.6) 

which was used in deriving equation (fTTl) . We first calculate R^ by using equation (lA.ip , 

R^ = ig - 2sin^9S + [2(1 - cos t9) + sin^ 

— 2 sin-t? (l — cos 79)8^ + (1 -cos^9)^S^ (A.7) 

Noting that tr (S^) = and tr (S^) = 2 and using equation flA.4l) . we obtain the trace 
of R2, 

tr (R2) = 4 cos^ ?9 - 1 = tr (R)^ - 2 tr (R). (A.8) 
This proves the relation in equation (]A.6p . 



Appendix B. Proof of equation (1261) 

In order to prove equation fl26l) . we recall the definition of Cki as in equation fl25|) . 
Substituting equation ( l22l) into equation ( 125|) . we get 

1 



C, 



kl 



dsL (a-^Rfca) (a^R^a) — J dadh (a'^R^a) (b'^R^b) 

T 1 

da (a (g) a)^ (Rfc O R/) (a O a) tr (R^) tr (R;) 

9 



(B.l) 



Let us define a couple of quantities, 

Ai = tr (RfeRf + RfcR;), 
A2 = tr(Rfc) tr(R;), 

^3= tr(Rfc)+tr(RO. (B.2) 

Then we rewrite the first term in equation fIB.ll) by using Schur's lemma, as in equation 
ffT6|l . so that 



da (a ® a) (R^ (g) R;) (a ® a) 



1 
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tr (Rfc) tr (R) + tr (RfcRH + tr (RkRi 



Then, Cki of equation (IB.ip is reduced to 
3Ai - 2A2 



kl 



6^ X 5 



- {A, + A2) . (B.3) 



(B.4) 



Using equation (lA.ip . Ai is explicitly calculated: 



cos-i^fc cost?/ + (cos-i^fc + COS??;) + n/| (1 — cosi^fc) (1 — cost? 



(B.5) 



Noting the last (third) term in equation ( IB.SP is semi-positive, is upper bounded by 

^l,max = 4 COS COS + 2 = - A3 + 3. (B.6) 
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It reaches the upper bound max when two rotation axes and are parallel or 
anti-parallel, i.e. |n^n;| = 1. On the other hand, Ai is lower bounded by 

^l,min = 2 (cos l^k COS l^i + COS l^k + cosi^i) = ^ {A2 + A3 - 3) . (B.7) 

The lower bound Ai^min is reached when and are ortogonal to each other or 
n^n, = 0. 

By substituting equation (IB.6P into equation (IB.4p . the upper bound of Cki is given 

by 

_ A2 - 3^3 + 9 _ [3-tr(Rfc)][3-tr(R;)] 
6^ X 5 6^ X 5 

Reminding of Aiq = [3 — tr (R)]/6-\/5 in equation (fT7|) . 

Using equation (IB.Tp . similarly, the lower bound of Cki is given by 



In case of = /, jn^n^l = 1 and thus 



Cki = ^Iq,u- (B.IO) 
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